Exact Spin and Pseudo-Spin Symmetric Solutions of the 



Dirac-Kratzer Problem with a tensor potential via Laplace 



Exact bound state solutions of the Dirac equation for the Kratzer 
potential in the presence of a tensor potential are studied by using the 
Laplace transform approach for the cases of spin- and pseudo-spin sym- 
metry. The energy spectra is obtained in the closed form for the rela- 
tivistic as well as non-relativistic cases including the Coulomb potential. 
It is seen that our analytical results are in agrement with the ones given 
in literature. The numerical results are also given in a table for different 
parameter values. 

1 Introduction 

The analytical solvable problems within the framework of quantum mechanics 
are very restricted. In this manner, it could be interesting to find the exact an- 
alytical solutions of the Kratzer potential in the relativistic domain. We intend 
to study the radial Dirac equation for the potential having the form 



where D > and a > 0. The detailed analysis of the bound states of this po- 
tential including a third parameter q in the non-relativistic domain are placed 
in literature [HH3]] where used the WKB approximation and the semiclassical 
quantization procedure with the help of the Langer transformation. The poten- 
tial form including the parameter q provides the extension of the spectra of the 
usual Kratzer potential 0ED . 

We find the exact bound state solutions and the corresponding upper- and 
lower-spinor of the radial Dirac equation in the presence of a tensor interac- 
tion within the framework of the Laplace transformation (LT) scheme [|6lf9l1 
for the cases of spin symmetry S(r) = +V(r), and pseudo-spin symmetry, 
S(r) = — V(r). The studying the solutions of the Dirac equation for the spin 
and pseudo-spin symmetric cases has been received a great attention in litera- 
ture IHPHI2H. 

The organization of this work is as follows. In Section 2, we find the energy 
spectrum and corresponding wave functions of the Kratzer potential II 13141 611 
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for the spin and pseudo-spin symmetric cases which have been studied in de- 
tails for the above potential without the tensor interaction in Ref. [17]. We also 
study the energy spectra of the Coulomb potential and see that the result ob- 
tained for the Coulomb problem is the same with the one obtained in Ref. [ 18] . 
At the end, we give the bound state spectrum of the above potentials for the 
non-relativistic limit where the results are consistent with the ones obtained in 
the literature. 

2 Bound States 

Dirac equation for the scalar S(r) and vector V(r) potentials with a tensor 
interaction U(r) is (ft = c = 1) QT9J 

{a.p + p [m +S(r)] + V(r) - ijB.fI7(r) - E nK }W) = 0, (2) 

where p is the momentum operator, E nK is the relativistic energy of the particle 
with the rest mass m and a written in terms of Pauli matrices and /3 are 
4x4 Dirac matrices, respectively. The energy of the particle E nK depends on 
the radial quantum number n and spin-orbit quantum number k. Inserting 
the eigenfunctions in terms of the spherical harmonic functions Y K (6,(j)) and 




F{r)Y K {e,ct>) 



(3) 



gives the following coupled differential equations for the upper-spinor F(r) and 
lower-spinor G(r) 

(4a) 
(4b) 

Substituting Eq. (4a) into Eq. (4b), we obtain a second-order differential equa- 
tion 

d 2 k(k + 1) d 2k j-[S(r)-V(r)][j- + Z-U(r)] 

-[- + U(r)]U(r)- dl 



d 


K 




+ - 


dr 


r 


d 


K 


dr 


r 



F(r) 



dr 2 r 2 L dr r m + E nK + S(r) — V(r) 

~ [ m l - e Ik + 2 ( m oS(r) + E nK V(r)) + S 2 (r) - V 2 (r)] F(r) = . (5) 

We consider the case where the scalar and vector potentials are equal, 
S(r) = V(r), which corresponds to the spin symmetry of the Dirac equa- 
tion I110H12II . Inserting Eq. (1) into Eq. (5), taking the tensor interaction 
as —C/r with a constant C and using a new wave function as F(r) = ■Jr 0(r) 
gives 

d 2 d>(r) 1 dd>(r) 1 , „ „ „ „>. 
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where 

a 2 = K ( K + 1) + C(l + 2k + C) + 2Da 2 (m + E n J + - , (7a) 
a 2 = -4L>a(m + E n J, (7b) 

Defining a wave function <^>(r) = r A /(r) with A as a constant and inserting into 
Eq. (6) leads 

d 2 /(r) d/(r) , , , . 

- (2a! - 1)-^ - (a 2 + e 2 K r) / (r) = 0, (8) 

where we set A = — to obtain a finite solution when r — » oo. By using the LT 
defined as 11201 



/• 00 

tr 



dre- tr f(r), (9) 
o 



we obtain a first-order differential equation from Eq. (8) as 

(t 2 - 4c) ^JT + [ (2a a + 1) t + a 2 ] f(t) = . (10) 
It's solution is 

/(t)-(t + e nif )-^ +1) [^] - " . (11) 



- hk 



a 2 2aj+l 



We impose the following condition to obtain a single-valued wave function 



a\ 1 



2e nK . 2 



-(2a a + l) = n. (n = 0,l,2,3,...) (12) 



Taking into account this condition and applying a simple series expansion to 
Eq. (11) gives 

" f_lW n | 

W ~ S r V. (2e - )m ( f + e -)" (2ai+1)_m > ^ 
i (n — m)!m! 

m=0 

Using the inverse Laplace transformation [20] in Eq. (13) we obtain 

A (-l) m n! r(2ai+l) 

/(r)~r 2ai e- e «- r V -4 ^- — — ^-(2e n(C r) m , (14) 

£ J («-m)!m!r(2a 1 + l + m)' " K 

and 

F(r) - iVr ai+ 2 e - e »- r ^^-n,!^ + l,2e nK r), (15) 
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where N is a normalization constant and used the following identity for the 
hypergeometric functions II21II 



1 F 1 (-n,a,x) = V — 



^ (-l)Pn! r(a) 



^(n-p)!p!r(cj+p) 



(16) 



The relation between the Laguerre polynomials and confluent hypergeometric 
functions as L^(x) = n, 17 + 1> X ) [21] gives the wave functions 

as 



Fir)-- 



n!r(2a-,+l) ,1 ,,, . 

n + 2a x + 1 



r ai+ 5e- e "* r I (2ai) (2e nK r). 



(17) 



Inserting the parameters in Eq. (7) into Eq. (12), we obtain the energy spectra 
of the Kratzer potential in the presence of a Coulomb-type tensor interaction 



m 



1 + 2n + ^(1 + 2k) 2 + 4C(1 + 2k + C) + 8Da 2 (£ nK . + m ) 



16D 2 a 2 



1 + 2n + ^(1 + 2k) 2 + 4C(1 + 2k + C) + 8Da 2 (£ nK . + m ) 



•(18) 



+ 16D 2 a 2 



which is the same with the ones obtained in Ref. II17II for C = 0. In Table 1, 
we present the numerical results for the case where C to see the effect 
of the tensor interaction on the energy spectrum. It is seen that the energy 
eigenvalues decrease while the strength of the tensor interaction increases. It 
is interesting to obtain the solution for the Coulomb case, setting the second 
term in Eq. (1) to zero, where we set the parameters as Da = Ze 2 (Ze 2 is the 
charge of the nucleus) and k = I in Eq. (18). This gives in the absence of 
tensor potential (h = c = 1) 



8Z z e 



2„4 



(n + l + l) 2 + 4Z 2 e 4 



(19) 



It is worth to say that this result is the same with the one obtained for the 
well-known Dirac-Coulomb problem in Ref. [18]. 

Now we intend to give the results for the Kratzer potential and Coulomb 
potential in the non-relativistic limit. 

2.1 Non-Relativistic Kratzer Limit 

In order to obtain the energy eigenvalues of the Kratzer potential in the non- 
relativistic limit for C = 0, we set E nK — m c 2 ~ E n ^ and E nK + m c 2 ~ 2m c 2 . 
Using this assumptions and with the help of Eq. (7), Eq. (12) gives (h = c = 1) 



8D 2 a 2 m 



-■nl 



1 + 2n + \ si 1 + 4£(£ + 1) + 16Dm a 2 



(20) 



which is exactly same with the result obtained in Ref. [22] . 
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2.2 Non-Relativistic Coulomb Limit 



We expand Eq. (19) into a series in terms of the parameter Ze 2 in order to 
obtain the energy spectrum of the Coulomb potential in the non-relativistic 
limit in the absence of tensor interaction. So, we find the energy equation for 

Ze 2 < 1 



8Z 2 e 4 



+ 



32Z 4 e 8 



(n + i + 1) 2 {n + i + lf 



+ ... 



(21) 



which is in agreement with the results obtained in literature [18]. 

Finally, we briefly study the case where the scalar and vector potentials 
are equal magnitude but different sign, i.e., S(r) = — V(r) in which the Dirac 
equation has pseudo-spin symmetry 111 0441 211 . For convenience for the rest of 
computation, we give the differential equation satisfying by the lower spinor- 
component. Following the same steps, we obtain the following 



d 2 (£(r) l#(r) 1 



dr 1 



+ - 

r 



-{A\+A\r + e 2 nK r 2 ) ( t>{r) = 0, 



where 



(22) 



A 2 = k(k - 1) + C{2k + C - 1) - 2Da 2 (m - E nK ) + - , 



Al = 4Da(m - E nK ) ; 



2 2 



E 2 =e 2 , 

TIK UK 



(23a) 

(23b) 
(23c) 



The wave functions can be written as 

n!r(2A 1 + l) i 

r l 1 ?. 



G(r) 



n + 2A 1 + 1 



-<Wr(2A 1 ) 



(2e„ K r) ; 



(24) 



and the energy eigenvalue equation of the Kratzer potential with a tensor in- 
teraction is written as 



E n K — m 



1 + 2n + a/(1 - 2k) 2 + 4C(2k + C - 1) + 8Da 2 (E nK - m c 2 ) 



1 + 2n + y/(l - 2k) 2 + 4C(2k + C - 1) + 8Da 2 (E nK - m c 2 ) 



+ 16D 2 a 2 



which is the same with the ones obtained in Ref. [17] C = 0. Our numeri- 
cal results are given for the pseudo-spin symmetric case in Table 1 for C ^ 
and see that the energy eigenvalues increase while the strength of the tensor 
interaction increases. 
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3 Conclusion 

We have exactly solved the radial Dirac equation for the Kratzer potential in the 
presence of a tensor interaction term by using the Laplace transform approach. 
The energy eigenvalues and the corresponding upper- and lower-spinor com- 
ponents of the potential are computed for the cases of the spin and pseudo-spin 
symmetry. We have also obtained the energy eigenvalues of the Coulomb po- 
tential by choosing suitable parameter values. We have also studied the non- 
relativistic limit for the above potentials. 
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Table 1: Energy eigenvalues of the Kratzer potential in the presence of the 
tensor interaction for the parameter values as m = 5 fm -1 , D = 1.25 fm -1 , 
a = 0.35fm- 1 EZJ. 



spin symmetry 


n 


K 


C = 0.25 


C = 0.50 





-2 


4.63346 


4.57251 




-3 


4.79992 


4.76888 




-4 


4.88031 


4.86515 




-5 


4.92171 


4.91357 


1 


-2 


4.82214 


4.80240 




-3 


4.88546 


4.87236 




-4 


4.92329 


4.91564 




-5 


4.94578 


4.94114 


pseudo-spin symmetry 


n 


K 


C = 0.25 


C = 0.50 


1 


-1 


-4.46315 






-2 


-4.83578 


-4.80233 




-3 


-4.90668 


-4.89445 




-4 


-4.93862 


-4.93245 


2 


-1 


-4.75526 






-2 


-4.90136 


-4.88615 




-3 


-4.93728 


-4.93064 




-4 


-4.95581 


-4.95208 
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